Phys. Rev. C 



Separated Structure Functions for Exclusive K^h. and K^TP Electroproduction at 

5.5 GeV with CLAS 

D.S. Carman,34 K. Park,^^ B.A. Raue/^ K.P. Adhikari,^^ D. Adikarani,29 M. Aghasyan/^ M.J. Amaryan,^^ 
M.D. Anderson,^'' S. Anefalos Pereira/^ M. Anghinolfi,^^ H. Avakian,^'' H. Baghdasaryan,^^'^^ J. Ball/ 
N.A. Baltzell/ M. Battaglieri/^ V. Batourine,34 1. Bedlinskiy,22 A.S. Biselli," J. Bono/^ S. Boiarinov,34 
W.J. Briscoe/5 W.K. Brooks,^^ V.D. Burkert,^^ A. Celentano/^ S. Cliandavar,^^ G. Charles/ P.L. Cole/^ 
M. Contalbrigo/^ O. Cortes/^ y. Crede/^ A. D'Angelo/^'^^ N. Dashyan/^ R. De Vita/^ E. De Sanctis/^ 
A. Deur/'' C. Djalali/^ D. Doughty/ Z'' R. Dupre/^ A. El Alaoui/ L. El Fassi/ P. Eugenio/^ G. Fedotov/^ 
S. Fegan/*^ R. Fersch/ J. A. Fleniing/° A. Fradi/^ M.Y. Gabrielyan/^ N. Gevorgyan/^ K.L. Giovanetti/^ 
F.X. Girod/'' J.T. Goetz/^ W. Gohn/ R.W. Gothe/^ K.A. Griffioen/s B. Guegan/i M. Guidal/^ L. Guo/^ 
/4 K. Hafidi/ H. Hakobyan/5'39 C. Hanretty/^ N. Harrison/ D. Heddle/'^* K. Hicks/^ D. Ho/ 
M. Hohrop/6 Y. Ilieva/3 D.G. Ireland/^ B.S. Ishkhanov/^ E.L. Isupov/^ H.S. Jo/^ K. Joo/ D. Keller/^ 
M. Khandaker/^ P. Khetarpal/^ A. Kim/" W. Kim/^ A. Klein/^ F.J. Klein/ S. Koirala/^ A. Kubarovsky/o 
V. Kubarovsky/4 S.V. Kuleshov/^'22 N.D. Kvahine/^ S. Lewis/" K. Livingston/" H.Y. Lu/ I .J 
.D. MacGregor/" Y. Mao/^ D. Martinez/" M. Mayer/^ B. McKinnon/" M.D. Mestayer/"* C.A. Meyer/ 

T. Mineeva/ M. Mirazita/^ V. Mokeev/** /2 R.A. Montgomery/" H. Moutarde/ E. Munevar/"* C. 
Munoz Camacho/^ P. Nadel-Turonski/'' R. Nasseripour/3'i2 c.S. Nepah/^ S. Niccolai/^ G. Niculescu/^ 
L Niculescu/3 M. Osipenko/^ A.L Ostrovidov/^ L.L. Pappalardo/^ R. Paremuzyan/^ S. Park/^ E. Pasyuk/** 
E. Phelps/3 J.J. Phillips/" S. Pisano/* O. Pogorelko/2 S. Pozdniakov/2 J.W. Price/ S. Procureur/ 
Y. Prok/'34 D. Protopopescu/" A.J.R. Puckett/'* G. Ricco/'* D. Rimal/2 M. Ripani/^ G. Rosner/" 
P. Rossi/4 Y. Sabatie/ M.S. Saini/^ C. Salgado/^ N.A. Saylor/" D. Schott/^ R.A. Schumacher/ E. Seder/ 

H. Seraydaryan/9 Y.G. Sharabian/^ G.D. Smith/" D.L Sober/ D. Sokhan/i S.S. Stepanyan/" 
S. Stepanyan/4 P. Stoler/o LL Strakovsky/^ S. Strauch/^ M. Taiuti/'' W. Tang/^ C.E. Taylor/" Y. Tian/^ 
S. Tkachenko/'^ A. Trivedi/^ M. Ungaro/^ B .Vernarsky/ H. Voskanyan/^ E. Voutier/'^ N.K. Walford/ 
L.B. Weinstein/^ M.H. Wood/ N. Zachariou/^ L. Zana/" J. Zhang/"* Z.W. Zhao/^ L Zonta/° 

(CLAS Collaboration) 
^ Argonne National Laboratory, Argonne, Illinois 60439 
^ Arizona State University, Tempe, Arizona 85287 
^ California State University, Dominguez Hills, Carson, California 90747 
^ Canisius College, Buffalo, New York 14208 
^ Carnegie Mellon University, Pittsburgh, Pennsylvania 15213 
^ Catholic University of America, Washington, D. C. 20064 
CEA, Centre de Saclay, Irfu/Service de Physique Nucleaire, 91191 Cif-sur-Yvette, France 
* Christopher Newport University, Newport News, Virginia 23606 
^ University of Connecticut, Storrs, Connecticut 06269 
Edinburgh University, Edinburgh EH9 3JZ, United Kingdom 
Fairfield University, Fairfield, Connecticut 06824 
Florida International University, Miami, Florida 33199 
Florida State University, Tallahassee, Florida 32306 
^* Universita di Genova, I6I46 Genova, Italy 
The George Washington University, Washington, D.G. 20052 
*® Idaho State University, Pocatello, Idaho 83209 
^"^ INFN, Sezione di Ferrara, 44100 Ferrara, Italy 
^'^ INFN, Laboratori Nazionali di Frascati, 00044 Frascati, Italy 
INFN, Sezione di Genova, I6I46 Genova, Italy 
INFN, Sezione di Roma Tor Vergata, 00133 Rome, Italy 
Institut de Physique Nucleaire ORSAY, Orsay, France 
Institute of Theoretical and Experimental Physics, Moscow, 117259, Russia 
James Madison University, Harrisonburg, Virginia 22807 
^* Kyungpook National University, Daegu 702-701, Republic of Korea 
2^ LPSC, Universite Joseph Fourier, CNRS/IN2P3, INPG, Grenoble, France 
University of New Hampshire, Durham, New Hampshire 03824 
■^^ Norfolk State University, Norfolk, Virginia 23504 
Ohio University, Athens, Ohio 45701 



Old Dominion University, Norfolk, Virginia 23529 
'^'^ Rensselaer Polytechnic Institute, Troy, New York 12180 

Universita' di Roma Tor Vergata, 00133 Rome, Italy 
Skobeltsyn Nuclear Physics Institute, 119899 Moscow, Russia 



2 



University of South Carolina, Columbia, South Carolina 29208 
Thomas Jefferson National Accelerator Facility, Newport News, Virginia 23606 
Universidad Tecnica Federico Santa Maria, Casilla 110-V Valparaiso, Chile 
University of Clasgow, Clasgow C12 8QQ, United Kingdom 
University of Virginia, Charlottesville, Virginia 22901 
College of William and Mary, Williamsburg, Virginia 23187 and 
Yerevan Physics Institute, 375036 Yerevan, Armenia 
(Dated: December 7, 2012) 

We report measurements of the exclusive electroproduction of K'^A and K^TP final states from 
an unpolarized proton target using the CLAS detector at the Thomas Jefferson National Accelerator 
Facility. The separated structure functions au, olt, o^tt, and cr^r/ were extracted from the 
dependent differential cross sections acquired with a longitudinally polarized 5.499 GeV electron 
beam. The data span a broad range of momentum transfers from 1.4 to 3.9 GeV^, invariant 
energy W from threshold to 2.6 GeV, and nearly the full center-of-mass angular range of the kaon. 
The separated structure functions provide an unprecedented data sample, which in conjunction with 
other meson photo- and electroproduction data, will help to constrain the higher-level analyses being 
performed to search for missing baryon resonances. 
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I. INTRODUCTION 

A complete mapping of the nucleon excitation spec- 
trum is the key to a detailed understanding of the ef- 
fective degrees of freedom of the nucleon and its asso- 
ciated dynamics. The most comprehensive predictions 
of this spectrum have come from various implementa- 
tions of the constituent quark model incorporating bro- 
ken SU(6) symmetry Additional dynamical contribu- 
tions from gluonic excitations in the wavefunction may 
also play a central role [2|] and resonances may be dynam- 
ically generated through baryon- meson interactions Q. 
Quark model calculations of the nucleon spectrum have 
predicted more states than have been seen experimen- 
tally |4]. This has been termed the "missing" resonance 
problem, and the existence of these states is tied in di- 
rectly with the underlying degrees of freedom of the nu- 
cleon that govern hadronic production at moderate ener- 
gies Q. 

Ideally we should expect that the fundamental theory 
that governs the strong interaction. Quantum Chromo- 
dynamics (QCD), should provide a reliable prediction of 
the nucleon excitation spectrum. However, due to the 
non-perturbative nature of QCD at these energies, this 
expectation has not yet been fully realized. There has 
been notable recent progress in calculations of QCD on 
the lattice that has led to predictions of the nucleon exci- 
tation spectrum with dynamical quarks, albeit with un- 
physical pion masses @ . Calculations with improved ac- 
tions, larger volumes, and smaller quark masses continue 
to progress. 

In parallel, the development of coupled-channel mod- 
els, such as those developed by the groups at Bonn- 
Gatchina 0,i], Giessen Jiihch and EBAC 
have made significant progress toward deconvoluting the 
nucleon spectrum. These multi-channel partial wave 
analyses have employed partial wave fits from SAID [T2| 
based on ttN elastic data to determine the properties of 



most N* and A* resonances listed in the Particle Data 
Group (PDG) [l^. Further critical information on the 
decay modes was obtained by including the inelastic re- 
actions ttN rjN, KA, K'E, and ttttA^. 

Recently the data landscape has undergone significant 
change with the publication of a vast amount of precision 
data in the photoproduction sector from JLab, SPring- 
8, MAMI, Bonn, and GRAAL. Data sets spanning a 
broad angular and energy range for jp — p7r°, mr^ , prj, 
Ptt^tf", pTT+TT", pn^r], K^A, and i^T+E" have provided 
high precision differential cross sections and polarization 
observables. Furthermore, new observables with polar- 
ized beams on both polarized proton and neutron targets 
have recently been acquired at several facilities and will 
be published over the next several years. 

In the K^A and if+E" electroproduction sector, 
dramatic changes to the world's database occurred 
with the publications from the CLAS Collaboration. 
These include (i) beam-recoil transferred polarization for 
K+A [U] and for i^+A and K+Y." 15], (ii) separated 
structure functions ajj = ar + ^ctl, o'lt, and cttt for 
K^A and ii'+E*', as well as ctt and [l6ll. and (iii) 
polarized structure function ctlt' for K'^A [17|. 

This paper now adds to and extends this database 
with the largest data set ever acquired in these kine- 
matics for polarized electrons on an unpolarized proton 
target. This work includes measurements of the sepa- 
rated structure functions ajj, olt, cttt, and ctlt' for 
the K^A and K^Y,^ final states at a beam energy of 
5.499 GeV, spanning W from threshold to 2.6 GeV, 
from 1.4 to 3.9 GeV^, and nearly the full center-of-mass 
angular range of the kaon. The full set of differential 
cross sections da/dfl*j^ included in this work consists of 
480 (450) bins in Q^, W, and cos 9*^ for the K+A {K+YP) 
final state and 3840 (3600) data points in Q'^, W, cos 6*^, 
and $ for K+A {K+TP). 

The organization for this paper is as follows. In Sec- 
tionini the different theoretical models that are compared 
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against the data are briefly described. In Section IIIIl 
the relevant formalism for the expression of the electro- 
production cross sections and separated structure func- 
tions is introduced. Section IIVI details the experimental 
setup and describes all analysis cuts and corrections to 
the data. Section |V] details the sources of systematic un- 
certainty on the measured cross sections and separated 
structure functions, which are presented in Section IVII 
along with a series of Legendre polynomial fits to the 
structure function data. Finally, we present a summary 
of this work and our conclusions in Section I VIII 



II. THEORETICAL MODELS 

To date the PDG lists only four N* states, 
iV(1650)l/2-, Af(1710)l/2+, Af(1720)3/2+, and 
7V(1900)3/2+, with known couphngs to KA and no iV* 
states are listed that couple to KT, [l^; only a single 
A* state, A(1920)3/2+, is listed with coupling strength 
to i^E. The branching ratios to KY provided for these 
states are typically less than 10% with uncertainties of 
the size of the measured coupling. While the relevance 
of this core set of N* states in the 7^*)^ K+A 
reaction has long been considered a well-established 
fact, this set of states falls short of reproducing the 
experimental results below W—2 GeV. Furthermore, 
recent analyses [l^ [l^ have called the importance of 
the iV(1710)l/2+ state into question. 

Beyond the core set of N* states, the PDG lists 
the iV(1900)3/2+ state as the sole established N* near 
1900 MeV. However, with a 500-MeV width quoted by 
some measurements, it is unlikely that this state by itself 
could explain the K^A cross sections below W=2 GeV, 
unless its parameters are significantly different than those 
given by the PDG. Recent analyses [13, [2l| have shown 
this state to be necessary to describe the CLAS beam- 
recoil polarization data ^2^. Note that the iV(1900)3/2+ 
state is predicted by symmetric quark models and its 
existence is not expected in diquark models. In the re- 
cent fits of "fp — if+S" data, all N* resonances found 
to be necessary to fit the K~^A data have been included. 
However, the existing database is smaller than the 

K~^A database, with significantly larger statistical uncer- 
tainties. 

A recent development in understanding the N* spec- 
trum was provided by the Bonn-Gatchina coupled- 
channel partial wave analysis of the hadronic ttN chan- 
nels and the photoproduced jp channels 0. This work 
presents an up-to-date listing of pole parameters and 
branching fractions for all N* and A* states up to 
^2 GeV with uncertainties at the level of a few percent. 
That analysis provided a list of (i) six N* states with cou- 
pling to KA, A^(1650)l/2-, iV(1710)l/2+, iV(1875)3/2-, 
iV(1880)l/2+, A^(1895)l/2-, iV(1900)3/2+, (ii) five N* 
states with coupling to KE, iV(1875)3/2-, iV(1880)l/2+, 
iV(1895)l/2-, Af(1900)3/2+, iV(2060)5/2-, and (iii) 
four A* states with coupling to iCE, A(1900)l/2-, 



A(1910)l/2+, A(1920)3/2+, A(1950)7/2+. For more on 
this list of states that couple to KA and ifE, see Ref. ^3] . 

The findings of Ref. are based on a significant 
amount of precision experimental data and the sophis- 
ticated coupled-channel fitting algorithms. However, 
in general, the issue of how to extract nucleon reso- 
nance content from open strangeness reactions is a long- 
standing question. Various analyses have led to very dif- 
ferent conclusions concerning the set of resonances that 
contribute (e.g. compare results from Refs. [2lj, and 
[2^ . as well as the statements made regarding the reso- 
nant set from Ref. Q)- Furthermore, lack of sufficient 
experimental information, incomplete kinematic cover- 
age, and underestimated systematics are still responsible 
for inconsistencies among the different models that fit 
the data to extract the contributing resonances and their 
properties Q. 

The indeterminacy for the open strangeness channels 
is in contrast to the pionic channels, where the contribut- 
ing resonances can be more reliably identified by means 
of a partial wave analysis for < 2 GeV. In open 
strangeness channels, this technique is less powerful as 
the non-resonant background contributions are a much 
larger fraction of the overall response. Several groups 
have stressed that the importance of the background con- 
tributions calls for a framework that accounts for both 
the resonant and non-resonant processes and that pro- 
vides for a means to constrain both of these classes of 
reaction mechanisms independently (26l . [27j . 

While there have been a number of publications of pre- 
cision cross sections and spin observables for both the 
photo- and electroproduction reactions, the vast major- 
ity of the theoretical effort has focused on fitting just the 
photoproduction data. Although KY photoproduction 
is easier to treat theoretically than KY electroproduc- 
tion, and is thus more amenable to a detailed quantitative 
analysis, the electroproduction reaction is potentially a 
much richer source of information concerning hadronic 
and electromagnetic interactions. The electroproduction 
observables have been shown to yield important com- 
plementary insights 26]. Some of the most important 
aspects of electroproduction include: 

• The data are sensitive to the internal structure 
of baryon resonances through the dependence 
of the electromagnetic form factors of the inter- 
mediate hadronic resonances associated with the 
strangeness production mechanism 3. 

• The structure functions are particularly powerful to 
gain control over the parameterization of the back- 
ground diagrams |28| . 

• Studies of finite processes are sensitive to both 
transverse and longitudinal virtual photon cou- 
plings, in contrast to the purely transverse response 
probed in the photoproduction reactions. 

• The longitudinal/transverse interference structure 
functions provide signatures of interfering partial 
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wave strengths that are often dramatic and have 
been shown to be useful for differentiating between 
models of the production amplitudes 



le 



13,^ 



• The beam-rccoil transferred polarizations in the 
fC+A and T^+E" reactions, as well as the recoil po- 
larization in the K reaction, have been shown to 
provide important new constraints to models that 
describe well the photoproduction data [3, HB] ■ 

At the medium energies of this work, perturbative 
QCD is not yet capable of providing predictions of dif- 
ferential cross sections. To understand the underlying 
physics, effective models must be employed that repre- 
sent approximations to QCD. Ultimately, it will be most 
appropriate to compare the electroproduction measure- 
ments against the results of a full coupled-channel partial 
wave analysis that is constrained by fits to the available 
data. Although output from such models is expected in 
the electroproduction sector in the future [sil |32| , as of 
now, these data have not yet been included in the fits. 
Thus comparisons of the electroproduction observables to 
single-channel models currently represent the best option 
to gain insight into the electroproduction realm. 

This analysis highlights three different theoretical 
model approaches. The first is a traditional hadrody- 
namic model and the second is based on K and K* 
Regge trajectory exchange. The third model, a hybrid 
Regge plus resonance approach, amounts to a cross be- 
tween the first two model types. Comparison of the dif- 
ferent model predictions to the data can be used to pro- 
vide indirect support for the existence of the different 
baryonic resonances and their branching ratios into the 
strange channels, as well as to improve constraints on 
the phenomenology of the different strangeness produc- 
tion reactions. The following subsections provide a brief 
description of the models included in this work. 




(a) 

s-channel 




(b) 
t-channel 



(c) 

u-channel 



FIG. 1: Tree-level diagrams contributing to the KY reactions: 
(a) s-channel exchanges, (b) t-channel exchanges, and (c) u- 
channel exchanges. 



Maxwell et al. [27|, |33|, [SJ] have developed a tree- 
level effective Lagrangian model (referred to as MX) 
for — >■ K'^A that incorporates the well-established 
s-channel resonances up to 2.2 GeV with spins up to 
5/2. The model also includes four A ?^-channel states, 
A(1405)l/2-, A(1670)l/2-, A(1820)5/2+, A(1890)3/2+, 
four E M-channel states, E(1385)3/2+, i;(1775)5/2-, 
S(1915)5/2+, E(1940)3/2-, and the K*{892) and 
ii'i(1270) f-channel resonances. 

The model was initially developed and fit to the avail- 
able 7p photoproduction data up to W—2.3 GeV 13411 . 
The most recent published version of the model [23| 
included fits to the available if +A separated structure 
function data from CLAS [l^. An extension of this 
model that also includes fits to the available CLAS K~^A 
(Jlt' data has been made available for this work as well. 
Overall the fits yield reasonable representations of both 
the photo- and electroproduction data. However, when 
compared to the results of a fit to the photoproduction 
data alone, the combined 7p and 7*p fit yields signifi- 
cantly different coupling parameters for an equally good 
overall fit to the data. This indicates that the pho- 
toproduction data alone are not adequate to uniquely 
constrain effective Lagrangian models of electromagnetic 
strangeness production. 



A. Hadrodynamic Model 

Hadrodynamic models provide a description of the re- 
action based on contributions from tree-level Born and 
extended Born terms in the s, t, and u reaction channels 
(see Fig. [T]) . The Born diagrams include the exchange of 
the proton, kaon, and ground-state hyperons, while the 
extended Born diagrams include the exchange of the asso- 
ciated excited states. This description of the interaction, 
which involves only first-order terms, is sensible as the in- 
cident and outgoing electrons interact rather weakly with 
the hadrons. A complete description of the physics pro- 
cesses requires taking into account all possible channels 
that could couple to the initial and final states, but the 
advantages of the tree-level approach are to limit com- 
plexity and to identify the dominant trends. The draw- 
back in this class of models is that very different conclu- 
sions about the strengths of the contributing diagrams 
may be reached depending on which set of resonances a 
given model includes. 



B. Regge Model 

Our KY electroproduction data are also compared to 
the Regge model from Guidal, Laget, and Vanderhaeghen 
[35I I (referred to as GLV). This calculation includes no 
baryon resonance terms at all. Instead, it is based only 
on gauge-invariant t-channel K and K* Regge-trajectory 
exchange. It therefore provides a complementary basis 
for studying the underlying dynamics of strangeness pro- 
duction. It is important to note that the Regge approach 
has far fewer parameters compared to the hadrodynamic 
models. These include the K and K* form factors and 
the coupling constants gxYN and gx'YN- 

The GLV model was fit to higher-energy photoproduc- 
tion data where there is little doubt of the dominance 
of kaon exchanges, and extrapolated down to JLab ener- 
gies. An important feature of this model is the way gauge 
invariance is achieved for the kaonic t-channel exchanges 
by Reggeizing the s-channel nucleon pole contribution in 
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the same manner as the t-channel diagrams. No counter 
terms need to be introduced to restore gauge invariance 
as is done in the hadrodynamic approach. 

The GLV Regge model reasonably accounts for the 
strength in the CLAS iiT+A differential cross sections 
and separated structure functions [l^ . Although the rea- 
sonable performance of a pure Regge description in this 
channel suggests a i-channel dominated process, there 
are obvious discrepancies between the Regge predictions 
and the data, indicative of s-channel strength. In the 
K+YP channel, the same Regge description significantly 
underpredicts the differential cross sections and sepa- 
rated structure functions [ll]. The fact that the Regge 
model fares poorly when compared to the K^Y!^ data is 
indicative that this process has a much larger s-channel 
content compared to K production. 



C. Regge Plus Resonance Model 

The final model included in this work was developed 
by the Ghent group [11], and is based on a tree-level ef- 
fective field model for i^+A and K^YP photoproduction 
from the proton. It differs from traditional isobar ap- 
proaches in its description of the non-resonant diagrams, 
which involve the exchange of K and K* Regge trajec- 
tories. A selection of s-channel resonances is then added 
to this background. This "Regge plus resonance" model 
(referred to as RPR) has the advantage that the back- 
ground diagrams contain only a few parameters that are 
tightly constrained by high-energy data. Furthermore, 
the use of Regge propagators eliminates the need to intro- 
duce strong form factors in the background terms, thus 
avoiding the gauge-invariance issues associated with tra- 
ditional effective Lagrangian models. 

In addition to the kaonic trajectories to model the 
t-channel background, the RPR model includes the 
same s-channel resonances as for the MX model be- 
low 2 GeV. The model does include several missing 
N* states at 1.9 GeV, A^(1900)3/2-, 7V(1900)3/2-, and 
A^(1900)l/2+. The separated structure functions [l6l . 
[itI and beam-recoil transferred polarization data from 
CLAS [l^ were compared to model variants with ei- 
ther a iV(1900)3/2- or a 7V(1900)l/2+ state at 1.9 GeV. 
Only the iV(1900)3/2^ state assumption could be rec- 
onciled with the data, whereas the A^(1900)l/2"'" option 
could clearly be rejected. In the K'^YP channel, four A* 
states, A(1700)3/2-, A(1900)l/2-, A(1910)l/2+, and 
A(1920)3/2+, have been included. 

In a new version of the RPR model (referred to as 
RPR-2011) [1^, several changes relative to the previous 
model version (referred to as RPR-2007) [1^ are note- 
worthy. The main difference is the implementation of an 
unbiased model selection methodology based on Bayesian 
inference. This inference is used as a quantitative mea- 
sure of whether the inclusion of a given set of N* states is 
justified by the data. Additionally, in this version of the 
model, the exchange of spin-3/2 resonances is described 



within a consistent interaction theory and the model has 
been extended to include the exchange of spin 5/2 reso- 
nances. 

The Regge background amplitude of RPR-2007 is con- 
strained by spectra above the resonance region {W > 
3 GeV) at forward angles (cosflj^ > 0.35). By extrap- 
olating the resulting amplitude to smaller W, one gets 
a parameter free background for the resonance region. 
The s-channel resonances are coherently added to the 
background amplitude. RPR-2007 describes the data 
for forward-angle photo- and electroproduction of K~^A 
and The resonance parameters of the RPR-2007 

model are constrained to the cos^J^ > 0.35 data. The 
RPR-2011 model with the highest evidence has nine well- 
established N* states and the "missing" states at 1.9 GeV 
with quantum numbers 7V(1900)3/2- and 7V(1900)l/2+, 
and has been fit to photoproduction data over the full 
center-of-mass (cm.) angular range. Neither ver- 
sion of the model has been constrained by fits to any of 
the electroproduction data. 



III. FORMALISM 

In kaon electroproduction a beam of electrons with 
four-momentum pe = (-Eg ,Pe) is incident upon a fixed 
proton target of mass Mp, and the outgoing scattered 
electron with momentum = {E(,i , pe' ) a-nd kaon with 
momentum pK = (Ek,Pk) are measured. The cross 
section for the exclusive K'^Y final state is then dif- 
ferential in the scattered electron momentum and kaon 
direction. Under the assumption of single-photon ex- 
change, where the virtual photon has four-momentum 
q = Pe — Pe' = (i^, q), this can be expressed as the prod- 
uct of an equivalent fiux of virtual photons and the "f*p 
cm. virtual photoabsorption cross section as: 



dEp'dVLp'cKl 



K 



, d^(Tv 

dfl* 



(1) 



where the virtual photon flux factor F depends upon 
only the electron scattering process. After integrating 
over the azimuthal angle of the scattered electron, the 
absorption cross section can be expressed in terms of 
the variables Q^, W, 9}^, and $, where = -Q^ 
is the squared four-momentum of the virtual photon, 

W = \J + 2Mpi' — is the total hadronic energy 

in the cm. frame, 9*j^ is the cm. kaon angle relative 
to the virtual photon direction, and $ is the angle be- 
tween the leptonic and hadronic production planes. A 
schematic illustration of electron scattering off a proton 
target, producing a final state electron, K'^ , and hyperon 
Y is shown in Fig. [2] 

Introducing the appropriate Jacobian, the form of the 
cross section can be rewritten as: 



dQ'^dWd^*j^ 



d^ay 

'dn*^' 



(2) 
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FIG. 2: Kinematics for K'^Y electroproduction defining the 
angles 9*jf and $. 



where 



a W W^-Ml 1 



(3) 



is the flux of virtual photons (using the definition from 
Ref. M), 



.= (l + 2-tan^- 



(4) 



is the polarization parameter of the virtual photon, and 
9f,i is the electron scattering angle in the laboratory 
frame. 

For the case of an unpolarized electron beam (helicity 
h=Qi) with no target or recoil polarizations, the virtual 
photon cross section can be written (using simplifying 
notation for the differential cross section) as: 



da 



{h = Q) = (To = <^u + £(Jtt cos 2$ 



K 



+ V'e(l + e)CTLT cos$, (5) 



where Ci are the structure functions that measure the 
response of the hadronic system and i — T, L, LT, and 
TT represents the transverse, longitudinal, and interfer- 
ence structure functions. The structure functions are, 
in general, functions of Q^, W, and 9*i^ only. In this 
work the unseparated structure function is defined as 
ajj — (7t + ecTL- 

In contrast to the case of real photons, where there is 
only the purely transverse response, virtual photons al- 
low longitudinal, transverse-transverse, and longitudinal- 
transverse interference terms to occur. Each of the struc- 
ture functions is related to the coupling of the hadronic 
current to different combinations of the transverse and 
longitudinal polarization of the virtual photon, ctt is 
the differential cross section contribution for unpolarized 
transverse virtual photons. In the limit — >■ 0, this 
term must approach the cross section for unpolarized real 



photons. (Ti is the differential cross section contribution 
for longitudinally polarized virtual photons, cttt and 
ctlt represent contributions to the cross section due to 
the interference of transversely polarized virtual photons 
and from transversely and longitudinally polarized vir- 
tual photons, respectively. 

For the case of a polarized electron beam with helicity 
h, the cross section form of Eq.© is modified to include 
an additional term: 



da 

dn^ 



ao + hy/e{l - e)aLT' sin 



(6) 



The electron beam polarization produces a fifth struc- 
ture function a^x' that is related to the beam helicity 
asymmetry via: 



Alt' = 



da 



e(l — f-)aLT' sin $ 



da 



(7) 



where the ± superscripts on correspond to the elec- 
tron helicity states of ft. = ±1. 

The polarized structure function ctlt' is intrinsically 
different from the structure functions of the unpolarized 
cross section. This term is generated by the imaginary 
part of terms involving the interference between longi- 
tudinal and transverse components of the hadronic and 
leptonic currents, in contrast to aLT, which is generated 
by the real part of the same interference. aLT' is non- 
vanishing only if the hadronic tensor is antisymmetric, 
which will occur in the presence of rescattering effects, 
interferences between multiple resonances, interferences 
between resonant and non-resonant processes, or even 
between non-resonant processes alone [sl]. ctlt' could 
be non-zero even when aLT is zero. When the reaction 
proceeds through a channel in which a single amplitude 
dominates, the longitudinal-transverse response will be 
real and ctlt' will vanish. Both aLT and aLT' are neces- 
sary to fully unravel the longitudinal-transverse response 
of the K^Y electroproduction reactions. 



IV. EXPERIMENT DESCRIPTION AND DATA 
ANALYSIS 

The measurement was carried out with the CEBAF 
Large Acceptance Spectrometer (CLAS) [s^ located in 
Hall B at JLab. The main magnetic field of CLAS is 
provided by six superconducting coils, which produce an 
approximately toroidal field in the azimuthal direction 
around the beam axis. The gaps between the cryostats 
arc instrumented with six identical detector packages. 
Each sector consists of drift chambers (DC) |4Ci for 
charged particle tracking, Cherenkov counters (CC) 41 1 
for electron identification, scintillator counters (SC) 42 1 
for charged particle identification, and electromagnetic 
calorimeters (EC) 43] for electron identification and de- 
tection of neutral particles. A 5-cm-long liquid-hydrogen 
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target was located 25 cm upstream of the nominal cen- 
ter of CLAS. The mam torus was operated at 60% of its 
maximum field value and had its polarity set such that 
negatively charged particles were bent toward the elec- 
tron beam line. A totally absorbing Faraday cup located 
at the end of the beam line was used to determine the 
integrated beam charge passing through the target. 

The efficiency of detection and reconstruction for sta- 
ble charged particles in the fiducial regions of CLAS is 
greater than 95%. The sohd angle coverage of CLAS is 
approximately Stt sr. The polar angle coverage for elec- 
trons ranges from 8° to 45°, while for hadrons it is from 
8° to 140°, with an angular resolution of 59, Scj) of better 
than 2 mr. The CLAS detector was designed to track 
particles having momenta greater than roughly 200 MeV 
with a resolution 6p/p of about 1%. 

The data in this paper were collected as part of the 
CLAS elf running period in 2003. The incident electron 
beam energy was 5.499 GeV. The live-time corrected in- 
tegrated luminosity of this data set is 10.6 fb~^. The 
data set contains 3.64 x 10^ e'K~^A events and 1.56 x 10^ 
e'i4r+S° events in the analysis bins included in this work. 

The data were taken at an average electron beam cur- 
rent of 7 nA at a luminosity of about 10^'' cm~^s~^. The 
event readout was triggered by a coincidence between 
a CC hit and an EC hit in a single sector, generating 
an event rate of ~2 kHz. The electron beam was lon- 
gitudinally polarized with polarization determined by a 
coincidence M0ller polarimeter. The average beam po- 
larization was about 75%. 

This analysis sought to measure the differential cross 
sections for the electroproduction reactions ep e'K^A 
and ep e'i<r+E° in bins of Q^, W, cos0|f, and 
Exploiting the $ dependence of the differential cross sec- 
tions (To as given by Eq.(l5]), a $ fit in each bin of Q^, W, 
and cos 0^ provides the separated structure functions au , 
(Tlt, and aTT- Finally, a $ fit to the beam spin asymme- 
try as given by Eq.d?]) in each bin of Q^, W, and cos^J^- 
gives access to the polarized structure function a^x'- 



A. Differential Cross Section Determination 

The bin-centered differential cross section for each hy- 
peron final state in each kinematic bin i was computed 
using the form: 



K 



1 



1 



(AQ2Avt/Acos0^A$) 
Rj ■ Ni ■ BCi 



(8) 



where F^, is the virtual photon flux factor computed ac- 
cording to Eq.(|31) for each bin at the bin-averaged mean 
of the bin and AQ^ AVFA cos 6'^A$ is the volume of each 
analysis bin computed using the bin sizes listed in Sec- 
tion |IVB] (the bin sizes are corrected for kinematic limits 
in the threshold W bins). Ri is the radiative correction 



factor, Ni is the background-subtracted ivT+A and K^TP 
yield in each bin, BCi is the factor that evolves the mea- 
sured bin-averaged differential cross section over each bin 
to a specific kinematic point within the Q^, W , cos0|f , $ 
bin, and rji accounts for the detector geometrical accep- 
tance and efficiency corrections. is the live-time cor- 
rected incident electron flux summed over all data runs 
included in this analysis determined from the Faraday 
Cup charge. For this experiment, the data acquisition 
live time ranged between 80 and 85%. The incident elec- 
tron flux was measured to be Nq = 9.807 x lO^^. Fi- 
nally, NApt/Au] represents the target number density, 
where Na is Avogadro's number, p=0. 07151 g/cm^ is 
the target density, t—b.Q cm is the target length, and 
yl^,=l. 00794 g/mol is the atomic weight of the target. 

The statistical uncertainty on the cross section in each 
bin i includes contributions from the statistical uncer- 
tainty on the hyperon yield and the acceptance function 
and is given by: 



(5(7,; = 



-I 1/2 



(9) 



B. Particle Identification and Event Selection 

The 7*p K+A and j*p — s> X+E" reaction channels 
were identified by detecting a scattered electron in co- 
incidence with a and then using the missing mass 
technique to identify the hyperons. Event reconstruction 
required the identification of both a final state electron 
and K'^ candidate within the well-understood fiducial re- 
gions of the detector. Details on the algorithms employed 
to minimize the particle misidentification at this stage 
are included in Ref. [ISi]. Before computing the missing 
mass spectrum, vertex cuts were employed to ensure that 
the particles originated from the target. In addition, cor- 
rections to the electron and kaon momenta were devised 
to account for reconstruction inaccuracies that arose due 
to to relative misalignments of the drift chambers in the 
CLAS magnetic field, as well as from uncertainties in the 
magnetic field map employed during charged track re- 
constructions. These corrections were typically less than 
1%. 

The algorithm used for hadron identification relied on 
comparing the measured velocity (3 = v/c for the track 
candidate to that expected for an assumed 7r+, , and 
p track. The assumption that resulted in the minimum 
A/3 — P — /3™j^ p was used to identify the species of the 
track. Fig. [3] shows A/3 versus momentum for the 
track assumption. For the data included here, the kaon 
momentum range was between 0.35 GeV (software cut) 
and w 4.5 GeV (kinematic limit), with a typical flight 
path of 5.5 m. The measured mass resolution was pri- 
marily due to the reconstructed time-of-flight resolution, 
which was «100 ps (a) on average; it also included con- 
tributions from the momentum and path length uncer- 
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tainties of CLAS. Fig. [3] shows that unambiguous sep- 
aration of K'^ tracks at the 2(T level is possible up to 
about 2 GeV. For higher momenta, the background due 
to particle misidentification increases. Detailed back- 
ground subtractions are necessary to determine the final 
event yields. 




FIG. 3: (Color online) A/3 vs. momentum (GeV) for the 
assumption that the reconstructed positively charged particle 
was a kaon. The K'^ band lies along A/3 — 0. 

Fig.|3]shows the e' missing mass {MM{e'K'^)) dis- 
tribution for the final event sample after all cuts have 
been made. This distribution contains a background con- 
tinuum beneath the hyperons that arises due to multi- 
particle final states where the candidate K'^ results from 
a misidentified pion or proton. 
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FIG. 4: (Color online) Distribution of MM{e'K+) (GeV) 
showing the A, TP, and several low- lying excited hyperon 
states. These data for the final event sample highlight the 
hyperon yields relative to the underlying background. The fit 
of the A peak shows the average mass resolution of a=15 MeV 
for this analysis. 

The data were binned in a four-dimensional space of 
the kinematic variables Q^, W, cosO}^, and $. The bin 
definitions used in this analysis are listed in Table HI Fig.[S] 



shows the kinematic extent of the data in terms of 
versus W and $ versus cos^J^-. These plots are overlaid 
with a grid indicating the bins in this analysis. The bin 
widths in W and $ were chosen to be uniform. Note that 
the maximum W bin at each was limited to where 
the hyperon yield fits were not dominated by systematic 
uncertainties. 
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FIG. 5: (Color online) Kinematic extent of the CLAS elf 
data set. (a). (GeV^) vs. W (GeV). (b). <1> (deg) vs. 
cosS^- The plots are overlaid with the binning choices in 
this analysis. 



C. Yield Extraction 

The three components of the MM{e'K^) spectra 
are the K events, the K^YP events, and the par- 
ticle misidentification background (dominated by pions 
misidentified as kaons). These individual contributions 
must be separated to extract the K'^K and K^Y}^ differ- 
ential cross sections in each analysis bin. 

The approach to separate the signal from the back- 
ground events employed a fitting process based on hy- 
peron template shapes and a polynomial to account for 
the particle misidentification background. The form for 
the spectrum fits was given by: 



MM ^ A- L 



template 



B ■ s 



template 



bck 1 



(10) 



where ^template and ^template are the simulated hyperon 
distributions with scaling factors A and B, respectively, 
and Pbck is a polynomial describing the background. 

The hyperon templates were derived from a GEANT- 
based Monte Carlo that included radiative processes 
and was matched to the detector resolution (see Sec- 
tion |TVDT]). The background contributions for this fit- 
ting were studied with a number of different assumptions 
(see discussion in Section |V]). Ultimately, a linear form 
for the background was chosen. The template fits to the 
missing mass spectra were carried out using a maximum 
log likelihood method appropriate for the statistical sam- 
ples of our data. Fig.[n]shows two sample fits to illustrate 
the typical fit quality to the data. 

The final yields in each kinematic bin were determined 
by taking the number of counts determined from the fits 
that fell within a mass window around the A (1.07 to 
1.15 GeV) and (1.17 to 1.22 GeV) peaks. Hyperon 
events in the tails of the distributions that fell outside of 
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g^: [1.4,2.2 GeV^] I W : [1.6,2.6 GeV] (20 50-MeV-wide bins) 
[2.2,3.0 GeV^] W : [1.6,2.4 GeV] (16 50-MeV-wide bins) 
[3.0,3.9 GeV^] | W : [1.6,2.2 GeV] (12 50-MeV-wide bins) 
cos 611, : [-0.9,-0.65], [-0.65,-0.4], [-0.4,-0.2], [-0.2,0.0], [0.0,0.2], 

[0.2,0.4], [0.4,0.6], [0.6,0.75], [0.75,0.9], [0.9,1.0] 

8 bins 45°-wide [-180°, 180°] 



TABLE I: Bin limits used for the KY cross sections and structure function analysis in this work. 




FIG. 6: (Color online) Sample template fits to the 
MM{e'K+) data (GeV) integrated over and $ for 
cos = 0.10 and IV=1.725 GeV (left) and 1.925 GeV (right) 
to illustrate the typical fit quality. The fit includes a A tem- 
plate, a template, and a polynomial background term. 

the mass windows were accounted for by the acceptance 
and radiative corrections. 

The number of A and S° hyperons in both the K~^A 
and iiT+E" mass windows relative to the total number of 
counts in the mass windows was found to be independent 
of and <I> in each bin of W and cos 9}^. Thus the 
final yields in each bin were determined by scaling the 
raw yields in the A'+A and K^Y.'^ mass windows by a 
background factor determined from fits in each bin of W 
and cos6't.. 



D. Acceptance and Efficiency Corrections 

1. Monte Carlo Acceptance Function 

Monte Carlo simulations were carried out for this anal- 
ysis for four distinct purposes. The first was to determine 
the detector acceptance in each bin, the second was as a 
cross check of the radiative correction factors, the third 
was to generate the hyperon templates for the spectrum 
fits, and the fourth was to determine the tracking effi- 
ciency corrections. 

For this analysis we employed two different event gen- 
erators for the exclusive if+A and X+E^ event samples. 
The first generator, FSGEN [11], generates ep e'K+Y 
events according to a phase space distribution with a 
t-slope scaled by a factor of e"''*. This generator did 
not include radiative effects. The nominal choice of the 
t-slope parameter of 6=1.0 GeV~^ was chosen to best 



match the cos 6*1- dependence of the data. The generated 
data were then weighted with ad hoc functions so that 
they matched well to the kinematic distributions of the 
data (see Fig.H]). 

The second generator, GENEV |45|, generates events 
for various meson production channels. It was modified 
for this analysis to include the i^+A and X+E^ chan- 
nels, reading in cross section tables for K^A and /C+E^ 
photoproduction based on the data of Refs. |;46j] and [13], 
respectively. It extrapolates to finite by introducing 
a virtual photon flux factor and electromagnetic form 
factors based on a simple dipole form. Radiative effects 
based on the formalism of Mo and Tsai [i^ are part of 
the generator as an option. Here too, the input distribu- 
tions of the model were weighted with ad hoc function so 
that they matched the data (see Fig.[7]). 




FIG. 7: (Color online) Comparison between selected A 
kinematic distributions {Q^ (GeV'^), W (GeV), cos^l, and 
"I" (deg)) of the data (black points with error bars) and the 
corresponding distributions generated from the FSGEN (red 
- light) and GENEV event generators (blue - dark). 

The Monte Carlo suite is based on a GEANT-3 pack- 
age [i^. The generated events were processed by this 
code based on the CLAS detector. The events were then 
subjected to additional smearing factors for the tracking 
and timing resolutions to match the average experimental 
resolutions. The analysis of the Monte Carlo data used 
the same code as was used to analyze the experimental 
data. Ultimately more than 1 billion Monte Carlo events 
were generated to determine the correction factors and 
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the associated systematic uncertainties, which are dis- 
cussed in Section IVl 

In order to relate the experimental yields to the cross 
sections, we require the detector acceptance to account 
for various effects, such as the geometric coverage of the 
detector, hardware and software inefficiencies, and reso- 
lution effects from the track reconstruction. The accep- 
tance is defined separately for the A'+A and K^Yf re- 
action channels as a function of the kinematic variables 



as: 



(11) 



where Nl'^'^ is the reconstructed number of events in each 
bin and Nf'^^ is the generated number of events in each 
bin. The FSGEN simulation was used to determine the 
acceptance function for the final analysis. Typical ac- 
ceptances for CLAS for the e' final state vary from 
wl% to 30%. Fig. [5] shows examples of this computed 
acceptance for the K^K final state as a function of $ 
and cos 6*}^ for one and W bin. 
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FIG. 8: Distribution of the computed K acceptance for 
CLAS as a function of cose^ and $ for the VF=1.925 GeV 
and (3^=1.8 GeV^ bin. The substructure in the acceptance 
is to due to the geometry of the active areas of the CLAS 
detector. The statistical error bars from the Monte Carlo are 
smaller than the symbol size on this plot. 



2. Efficiency Corrections 

For this analysis several standard CLAS efficiency cor- 
rections were applied to the yields on an event-by-event 
basis. The first correction accounted for the efficiency 
of the Cherenkov counter for registering electron tracks 
based on the number of detected photoelectrons in each 
sector in a fine grid of the 9 and angles of the electron 
at the face of each CC detector. The average CC effi- 
ciency within the electron geometric fiducial cuts for this 
analysis is 96%. 

The remaining efficiency corrections account for 
hadron tracking inefficiencies. The first correction ac- 
counts for the single track reconstruction efficiency in 
CLAS that is not 100% due to inefficient SC paddles 
and DC tracking regions. This efficiency function was 



assigned based on the relative ratio of data counts to 
Monte Carlo counts as a function of CLAS sector and 
SC paddle number. These corrections are at the level of 
about 10% on average. 

Another efficiency correction related to tracking is nec- 
essary for events in which two charged tracks of the same 
charge and similar momenta lie very close to each other. 
For such events the tracking algorithm may not success- 
fully identify two separate tracks. For this analysis, a 
correction was applied to the small fraction of events in 
which the K'^ and p from the decay of the A were in the 
same CLAS sector within 10° of each other in polar an- 
gle. This efficiency factor is necessary even for the e' 
analysis due the presence of the decay protons in the fi- 
nal state. The systematics associated with each of these 
efficiency corrections are discussed in Section IVBl 



E. Radiative Corrections 

Radiative effects must be considered when determining 
the 7*p K^Y cross sections. Radiative effects result 
in bin migration such that the measured and W are 
not the true and W to which the event should be 
properly associated. 

For this analysis, two different approaches to deter- 
mine these correction factors have been employed. The 
first uses the stand-alone program EXCLURAD [s^] and 
the second uses the event generator GENEV in com- 
bination with the CLAS Monte Carlo. The radiative cor- 
rection factor that multiplies the measured bin-averaged 
differential cross section in each bin is defined as the ratio 
of the computed bin-averaged cross section with radia- 
tion off to that with radiation on. More details on each 
program are included below. 



1. EXCLURAD 

EXCLURAD represents a covariant technique of can- 
cellation of the infrared divergence that leads to in- 
dependence of any parameter that splits the soft and 
hard regions of phase space of the radiated photons. 
It uses an integration technique that is exact over the 
bremsstrahlung photon phase space, and thus does not 
rely on the peaking approximation |5lj . This approach 
is an exact calculation in that it specifically accounts for 
the exclusive nature of the reactions as the detection of 
hadrons in the final state, in addition to the electron, 
reduces the phase space allowed for the final radiative 
photons. 

The program EXCLURAD was based on the mea- 
sured structure functions from this analysis for K and 
K'^YP . The structure functions uu, cflTi o'tt, and <tlt' 
were read into the program and the cross section ratio for 
each bin in Q^, W, cos 6*^, and $ was computed with ra- 
diation off to that with radiation on, giving the radiative 
correction factor Ri for that bin. 
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The trends of the correction (shown in Fig. 19]) are 
such that it has its largest value near threshold and then 
quickly falls off to a near constant average value with in- 
creasing W. Note that the radiative correction factors 
including the helicity-dependent structure function ctlt' 
for the two helicity states have no impact on the helicity 
asymmetry computation in Eq.Q and are not included 
in the analysis. 
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FIG. 9: (Color online) Radiative corrections factors for 
the K+A final state from EXCLURAD (blue - circles) and 
GENEV (red - triangles) for Q^=1.80 GeV^ as a function 
of cos0^ for representative values of W of 1.725, 1.925, and 
2.125 GeV. The GENEV calculations are independent of <1> 
and only one data point is shown at each value of cos 9^. 
The EXCLURAD calculations have a $ dependence that is 
symmetric about <1>=0. The different radiative correction val- 
ues for each <1> at a given cos6lJf are included on the plot. 



2. GENEV 

The event generator GENEV [i^ was introduced in 
Section IIVD II as it was used to compute the CLAS ac- 
ceptance function. This program also allows for radiative 
correction factors to be determined. It includes radiative 
effects based on the formalism for inclusive electron scat- 
tering from Ref. (48j and employs the peaking approx- 
imation [5l| in the computation. As GENEV is based 
on an evolution of the photoproduction cross sections, it 
does not have an explicit $ dependence and thus the Ri 
factors in Eq.® were determined in bins of Q^, W, and 
cos 9*j^. 



This model has several shortcomings. The first is that 
the phase space for the radiated photons is not properly 
computed as this is modified by the detected hadrons. 
Secondly, the model is based on only the longitudinal and 
transverse response and does not include the interference 
structure functions a^x or aTT- Finally, the approach 



relies on an unphysical parameter to split the hard and 
soft regions of the radiated photon phase space to cancel 
the infrared divergence. Due to the known limitations 
with this approach, it was used only to provide a qual- 
itative cross check to the EXCLURAD results and to 
explore the associated systematic uncertainties (see Sec- 
tion |VB|). Fig. m shows a comparison of the radiative cor- 
rection factors computed by GENEV to those computed 
from EXCLURAD. Apart from the region near thresh- 
old, the correspondence between the two approaches is 
within 10%. 



F. Bin Centering Corrections 

The goal of this analysis is to measure cross sections 
and separated structure functions for the K'^Y final 
states at specific kinematic points. However, the analysis 
proceeds from using finite bins in the relevant kinematic 
quantities Q^, W, cos 6'^, and $ (see Section HV Bp . 

The virtual photon flux factor F^ defined in Section Hill 
is computed for each bin using the bin-averaged values of 
and W. If the cross sections were computed at this 
point using Eq.® with the Bd terms set to unity, we 
would have completed a measurement of the bin-averaged 
cross sections that we could quote at the corresponding 
bin-averaged kinematic points. To quote the cross section 
at specific kinematic points of our choosing, namely, the 
geometric centers of the defined bins, we must evolve the 
cross sections from the bin-averaged kinematic points to 
the geometric bin centers. These evolution factors are 
the bin-centering correction factors BCi in Eq.®. The 
bin-centering corrections are then applied for each bin as: 



point 



, avg / daP°"^*" 



dn 



model 



(12) 

where BCi are the ratios of the bin-centered cross section 
to the bin-averaged cross section. 

Studies of the bin-averaged kinematic quantities versus 
the geometric bin-centered values show that there is no 
need for bin-centering corrections in W or cos0^. For 
this work the threshold W bin for if+A is quoted at 
1.630 GeV and for if+E" at 1.695 GeV. To determine 
the bin-centering factor BCi for each bin, we have fit the 
measured structure functions ajj for each W and cos 9}^ 
bin versus for both the K^A and if+E° final states. 
To bin center the data at specific points, we have used 
the following dipole evolution factor: 



BC, 



(1 



^ point 



/0.7) 



(l + QLg./O.T) 



(g2 in GeV). 



(13) 



The bin centering factors using this form were in the 
range from 0.95 to 1.05 across the full kinematic phase 
space. 
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G. Structure Function Extraction 

The differential cross sections computed using Eq.([S]) 
are the mean values within the finite size of the $ bins 
and therefore do not reflect the value at the bin center. 
Thus directly fitting these data with Eq. ^ to extract the 
structure functions ajj = ut + cctl, ctt, and ctlt would 
be inappropriate. Integrating Eq.([5]) over the finite bin 
size, A$ = $„ — <i>;, where and are the upper and 
lower limits of the bin, respectively, gives: 



ao = I {cru + ^CTTT cos 2$+ 



Ve(e+T)(TLT cos ^)d^ 

= -^A'^u^^ + ^CTTT (sin2'I>„ - sin 2$, 
A* 2 

y/e{e + l)aLT (sin - sin $;)). (14) 



ctq now represents the value of the measured bin-averaged 
cross section in a given $ bin and fitting the data with 
Eq. ([T4| yields the separated structure functions. The "e" 
pre-factors were evaluated at the bin center and divided 
out. Note that prior to the $ fits, the statistical uncer- 
tainty on each cross section point was combined linearly 
with that portion of the systematic uncertainty arising 
from the yield extraction procedures (see Section IV Al for 
details). 
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FIG. 10: Fits of the K'^A differential cross section (nb/sr) 
vs. $ (deg) for VK=1.725 GeV and Q^=1.8 GeV^ showing a 
sample of the fits to extract au, o^lt, and arr- 

In Fig. [10] we show a sample of the ^-dependent 
differential cross sections for the K^A final state at 
W=1.725 GeV for GeV^. The different shapes 

of the differential cross sections versus $ in each of our 
bins in Q^, W, and cos 6*^ refiect differences of the in- 
terference terms, a^T and axT, while the differences in 
scale reflect the differences in ajj- 
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FIG. 11: Normalized asymmetries for K^A vs. <& (deg) at for 
VF=1.725 GeV and g^=1.8 GeV^ showing a sample of the fits 
to extract ai^x' ■ 



The extraction of ctlt' in each bin of Q^, W, and 
cos^Jf requires knowledge of both the asymmetry A^t' 
and the unpolarized cross section ctq, which can be 
seen by rearranging Eq.(I7]) into a normalized asymme- 
try A'^'''"' as: 



A meas 



= cflt' sin $. 



(15) 



Alt' is determined by forming the asymmetry of the 
if+A and if+S" yields for the positive and negative 
beam helicity states {h = ±1) as: 



Alt' — 



1 

Pb 



N+ - N' 
N+ + N- 



(16) 



where Pb is the average longitudinal polarization of the 
electron beam. 

As with the cross sections, the measured asymmetries 
are the average values over the span of the given $ bins. 
Integrating Eq.([7]) over the size of the $ bin results in: 



A 



N 



A 7neas 



sin <I>A<1) 

cos <I>; — COS 



(17) 



To extract ctlt' , a sin $ fit was performed according to 
Eq. ([T7|) . where the kinematic e factor was calculated at 
the bin-centered values of and W for each bin. A sam- 
ple of these distributions is shown in Fig. [TT] for the K^A 
final state at 1^=1.725 GeV for Q^^l.S GeV^. Similar 
to the case for the unpolarized structure function extrac- 
tion discussed in Section IIV G| prior to the $ fits the 
statistical uncertainty on the helicity-gated yields was 
combined linearly with that portion of the systematic 
uncertainty arising from the yield extraction procedure 
(see Section IV Al for details). 

The statistical uncertainty on the data points in each 
bin i are a combination of the contributions from both 
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Alt' and ctq E^nd are given by: 

S{ALT"Jo)^ = \/{ALT'6ao)^ + {aoSALT')l (18) 

V. SYSTEMATIC UNCERTAINTIES 

To obtain a virtual photoabsorption cross section, we 
extract the yields for the A'+A and if +1]" reactions from 
the the missing-mass spectra for each of our bins in Q^, 
cos 9^, and The yields are corrected for the ac- 
ceptance function of CLAS including various efficiency 
factors, radiative effects, and bin-centering factors. Fi- 
nally, we divide by the virtual photon flux factor, the 
bin volume corrected for kinematic limits, and the beam- 
target luminosity to yield the cross section. Each of these 
procedures is subject to systematic uncertainty. We typ- 
ically estimate the size of the systematic uncertainties by 
repeating a procedure in a slightly different way, e.g. by 
varying a cut parameter within reasonable limits, by em- 
ploying an alternative algorithm, or by using a different 
model to extract a correction, and noting how the results 
change. 

In this section we describe our main sources of sys- 
tematics. The five categories of systematic uncertainty 
studied in this analysis include yield extraction, detector 
acceptance, radiative corrections, bin centering correc- 
tions, and scale uncertainties. Each of these categories is 
explained in more detail below. 

In assigning the associated systematic uncertainties, 
we have compared the differential cross sections and ex- 
tracted structure functions, ajj, a^T, ^-tt, and ct^t'i 
with the nominal cuts and the altered cuts. The frac- 
tional uncertainty for each bin i was calculated via: 



The relative difference in the results Sat is then used 
as a measure of the systematic uncertainty. In this anal- 
ysis we have carefully studied the kinematic dependence 
of the systematics and conclude that there is no evidence 
within a given bin of systematic variations with W, 
cos , or $. Tableimiists the categories, specific sources, 
and the assigned systematic uncertainties on our mea- 
surements. Overall the scale of the systematic uncertain- 
ties is at the level of about 10%. 



A. Yield Extraction 



detected tt"*" as a . We have concluded that all system- 
atic effects associated with the spectrum fitting get larger 
in direct proportion to the size of the statistical uncer- 
tainty. We estimated that the systematic uncertainty due 
to the yield extraction is roughly equal to 20% of the size 
of the statistical uncertainty in any given bin. We added 
these correlated uncertainties linearly with the statistical 
uncertainties on our extracted yields before performing 
the $ fits. 

The other sources of systematic uncertainty considered 
in this category are associated with the defined electron 
and hadron fiducial cuts and the cuts on the deposited 
energy in the calorimeter used to identify the candidate 
electron sample. Variations in the definitions of the fidu- 
cial cuts and the EC energy cuts over a broad range 
showed that the observables were stable for each cut type 
to within 5%. 



B. Detector Acceptance 

In the category of detector acceptance, the associated 
systematics include that due to the model dependence 
of the acceptance function, the stability of the tracking 
efficiency corrections, and the CC efficiency function. 

For this analysis both the FSGEN and GENEV physics 
models were used to generate the Monte Carlo events. 
Because of the finite bin sizes used in this analysis, it is 
necessary to study how the derived acceptance function 
based on the different event generators impacts the ex- 
tracted observables. For both models we determined the 
acceptance function and stepped through the full analy- 
sis chain to extract the observables. The systematics as- 
signed for the model dependence were in the range from 
about 4% to 9%. 

The approach to assign a systematic associated with 
the CLAS tracking efficiency corrections was to employ 
slightly different algorithms and then to step through the 
full analysis chain. The tracking efficiency gave stable 
results at the level of 5%. The systematic associated with 
the close track efficiency was stable in the range from 2 
to 5%. 

To study the systematic uncertainty associated with 
the CC efficiency function, we compared the measured 
observables with the nominal CC efficiency corrections 
to an analysis with the CC efficiency set to 100% for all 
events. The differences were within 1.5% for all observ- 
ables. 



The procedure to determine the K'^Y yields in each 
analysis bin employs hyperon templates derived from 
Monte Carlo simulations that have been tuned to match 
the data. The background fit function has been studied 
using two different approaches. The first uses a polyno- 
mial (either linear or quadratic) and the second uses the 
ep e'n^ X data sample purposefully misidentifying the 



C. Radiative Corrections 

Two very different approaches have been used to study 
the radiative corrections for the K'^ K and if+S^ elec- 
troproduction reactions. The first was the exclusive ap- 
proach based on the EXCLURAD program [s^l and the 
second was based on the inclusive approach based on the 
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dcr/dQ 1 au \o'lt: o'tt\ <^lt' 


1. Yield Extraction 




Signal fitting/binning effects 


1.20x Stat. err. 


Fiducial cuts 


0.4-2.6% 


0.4-2.6% 




0.7-4.4% 


Electron identification 


1.1% 


0.1% 


4.0% 


1.4% 


2. Detector Acceptance 




MC model dependence 


4.0-9.3% 


3.6-7.8% 


6.8% 


3.6-7.0% 


Tracking efficiencies 


5.3% 


5.3% 


5.5% 


5.3% 


Close track efficiencies 


2.8% 


1.6% 


4.7% 


2.6% 


CC efficiency function 


1.5% 


1.5% 


1.5% 


1.5% 


3. Radiative Corrections 


2.0% 


2.0% 


4.4% 


2.0% 


4. Bin Centering 


0.5% 


0.5% 


0.5% 


0.5% 


o. ocait u iicei LaiiiLich 












2.3% 


Ptioton flux factor 


3.0% 


3.0% 




3.0% 


Luminosity 


3.0% 


3.0% 




3.0% 


Total Q'i 


12.5% 


11.1% 


11.7% 


11.6% 


Total Q'i 


9.2% 


8.2% 


11.7% 


9.2% 


Total Q'i 


8.9% 


8.5% 


11.7% 


9.0% 



TABLE II: Categories and systematic uncertainty assignment for the observables reported in this work for our three Q^ points 

point are obtained by 



at Ql 



1.80, (32=2.60, and Q3=3.45 GeV . The total systematic uncertainty assignments for each Q 
adding the different contributions in quadrature. 



GENEV program |45| ■ Comparison of the extracted ra- 
diative corrections between EXCLURAD and GENEV 
were within about 8% of each other. However, due to 
the shortcomings of the GENEV model as discussed in 
Section Five 21 this comparison was only used as a cross 
check of the overall scale of the corrections. 

To assign a systematic uncertainty for the radiative 
corrections for this analysis, we compared the measured 
observables using the EXCLURAD approach but varying 
the energy range of integration of the radiated photon 
over a broad range. The corrections were stable in the 
range from 2 to 5%. 



D. Bin Centering Corrections 

To assign a systematic uncertainty to the bin center- 
ing corrections, the mass term in the dipole form (see 
Ea. p3p ) was varied over a broad range. The maximum 
variation seen in any of the extracted observables was 
0.5%. 



tematic assigned to ult'- This is given by: 



5a 



LT' 



meas 
LT' 



6P^ 



= Ictlt' 10.023, 



(20) 



where (5Pe=0.03 and Pf, = 0.754 is the average beam po- 
larization. Thus the assigned systematic for <tlt' due to 
the beam polarization uncertainty is 2.3%. 

The uncertainties in the average virtual photon flux 
factor across our phase space were estimated by propa- 
gating through the flux definition the uncertainties asso- 
ciated with W and that arise from the uncertainty in 
the reconstructed electron momentum and angles. The 
uncertainty in the flux factor was determined to be 3%. 
This scale-type uncertainty affects only the differential 
cross section and the structure functions ajj and ctlt'- 

We estimated uncertainties in the beam-target lumi- 
nosity based on the analysis of CLAS ep elastic scatter- 
ing cross sections from Ref. 52|. The overall systematic 



uncertainty of the Faraday Cup charge measurement has 
been assigned to be 3.0%. This scale-type uncertainty af- 
fects only the differential cross section and the structure 
functions au and ctlt'- 



E. Scale Uncertainties 



Cross Cliecks 



In the category of scale uncertainties, the associated 
systematics include that due to the beam-charge asym- 
metry and uncertainties in the beam polarization, the 
photon flux factor, and the luminosity. 

The estimated beam-charge asymmetry is at the level 
of a few times 10"'' and is thus entirely negligible. The 
uncertainty in the beam polarization affects only the sys- 



The nominal analysis for the K'^A and iiT+S^ differ- 
ential cross sections and separated structure functions 
required only the detection of the electron and if+ in 
the final state. In order to check the overall systematic 
assignment, the observables were also extracted when de- 
tecting an additional p. The detection of the proton 
from the A decay gives rise to an analysis sensitive to 
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the same systematic uncertainties as the nominal anal- 
ysis, and thus should yield consistent results. However, 
requiring the proton reduces the acceptance by roughly 
a factor of three, therefore this comparison can only be 
used as a cross check of the nominal analysis. 

The agreement between the cross sections extracted 
using the e'K^ and e'K^p final states is at the level 
of ±5-10% and independent of kinematics to within the 
statistical uncertainties. The differences are driven by 
the marginal statistics in some of the analysis bins for 
the e'K~^p analysis. These comparisons show that the 
assigned systematic uncertainties are reasonable. 

VI. RESULTS AND DISCUSSION 
A. Angular Dependence 

In Figs. [12] and [13] we show the extracted structure 
functions ajj, a^T, <^tt, and ctlt' versus cosflj, for the 
if+A final state. Figs . [HI and [T5l show the same plots for 
the K^YP final state. These plots are for our lowest 
point at 1.80 GeV^. The general conclusions that can be 
drawn from studying the angular dependence are simi- 
lar for the two higher points at 2.60 and 3.45 GeV^. 
However, the full set of our data is available in the CLAS 
physics database f53| . 

The following curves are overlaid on the data: 

• The hadrodynamic model of Maxwell et al. (MX) 
(red/dashed curves - thinner line type from 
Refs. [13, [131, thicker line type is an extension 
of that model including fits to glt' data from 
Ref. [13] )• Note that this model is only available 
for the JT+A final state and calculations go to a 
maximum W of 2.275 GeV. 

• The Regge model of Guidal et al. {GIN) [H 
(green/dotted). 

• The Regge plus resonance model of Ghent 
(RPR) [2^ (black/solid curves - RPR-2007 thin- 
ner line type, RPR-2011 thicker line type). For the 
K'^TP comparison, only the RPR-2007 version is 
presently available. 

A number of observations can be made independent of 
the model calculations: 

1. The production dynamics for K^h. and K^YP are 
quite different for W < 2 GeV. However, as W in- 
creases further, the production mechanisms become 
similar. This is to be expected as KY production 
is known to be dominated by i-channel exchanges 
at higher energies. 

2. The K~^A production dynamics are dominated by 
i-channel exchange over the full resonance region 
as indicated by the strong forward peaking of au 
in Figs. [T2] and [T3] However, given the mid-angle 



peaking of ajj for if+S" below 2 GeV, clearly s- 
channel contributions play a much more significant 
role for this final state. 

3. The forward peaking of au and (Tlt for K^A com- 
pared to K~^T}^ can be qualitatively explained by 
the effect of the longitudinal coupling of the vir- 
tual photons. We note that the two channels are 
of nearly equal strength at (5^=0 GeV^ [H, l47j . 
while here at Q^=1.80 GeV^ the K+A channel is 
stronger than the if+S^ channel at forward angles 
by a factor of 3 to 4. For transverse (real) photons, 
the i-channel mechanism at low t is dominated by 
vector K*'^ exchange, which relates directly to the 
magnitudes of the coupling constants gx'YN rela- 
tive to gKYN- As rises from zero, the photon 
can acquire a longitudinal polarization and the im- 
portance of pseudoscalar if+ exchange increases. 
Given that gj^j^j^ » g'kt.on |13,[5l], this effect in- 
creases the cross section for K^A relative to K^YP 
(this is consistent with the arguments presented in 
Ref. !l6]). This argument is consistent with our 
observation of a sizable glt for K'^ A and a glt 
consistent with zero for K'^YP . It should also be 
the case that since gK'Y.N ^ gx'SN, K* exchange 
should dominate the K^Y'^ channel. Because K* 
exchange must vanish at forward angles due to an- 
gular momentum conservation, the K^Y^ cross sec- 
tion should also decrease at forward angles [35j . 

4. For K^A, arr is consistent with zero up to 
about W—1.9 GeV then develops a strong for- 
ward peaking that abruptly changes sign at about 
W^=2.2 GeV. For K+Y°, utt peaks at mid-range 
angles up to W~2 GeV and then looks very similar 
to K'^A for higher W . This higher W response is 
well explained by the interference of the K and K* 
Regge trajectories. 

5. For iC+A, aLT' is relatively flat over the full an- 
gular range up to W=2 GeV and then develops a 
strong forward peaking for higher W very similar 
to the other interference structure functions. We 
also note that it is significantly reduced at this 
compared to the results at (5^=0.65 and 1.0 GeV'^ 
shown in Ref. [l7!|. ctlt' for K^Y'^ is consistent 
with zero over the full angular range. 

Comparing the data in Figs. [T^] to [15] to the differ- 
ent single-channel model calculations, it is apparent that 
none of the models is successful at fully describing all of 
the data. A few general remarks are in order: 

1. In general the models agree better with the K^A 
data than with the K^Y^ data. This likely arises, 
in part, due to the fact that better quality data for 
K^A is available than for if+E^. However, as the 
resonance content is stronger in if+E" compared to 
K^A for < 2 GeV given that the Regge predic- 
tions for A are in much closer agreement with 
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FIG. 12: (Color online) Structure functions au, ctlt, cttt, and ct^t/ (in nb/sr) for K'^A production vs. cosO'^ at 5.499 GeV for Q^^l.SO GeV^ and W from 1.630 
to 2.075 GeV. The error bars represent the statistical uncertainties only. The curves shown are from the model calculations of Maxwell et al. (MX) (red/dashed 
curves) [27l. [3^. |34|]. Guidal et al. (GLV) (blue/dot-dashed curves) =3^], and Ghent (RPR) (black/solid curves) [2^. See the text for detailed descriptions of the 
calculations and the corresponding references. 
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FIG. 13: (Color online) Structure functions au, (tlt, cftt, and ai^T' (in nb/sr) for K^A production vs. cos d^. at 5.499 GeV for Q^=1.80 GeV^ and W from 2.125 to 
2.575 GeV. The error bars represent the statistical uncertainties only. The curves are defined in the caption of Fig. 1121 




FIG. 14: (Color online) Structure functions au, olt, cttt, and cr^T' (in nb/sr) for K^TP production vs. cos6'J- at 5.499 GeV for Q^^l.SO GeV^ and W from 1.630 to 
2.075 GeV. The error bars represent the statistical uncertainties only. The curves are defined in the caption of Fig. 1121 




FIG. 15; (Color online) Structure functions au, olt, cttt, and cr^T' (in nb/sr) for K^TP production vs. cos6'J- at 5.499 GeV for Q^^l.SO GeV^ and W from 2.125 to 
2.575 GeV. The error bars represent the statistical uncertainties only. The curves are defined in the caption of Fig. 1121 
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the ajj measurements compared to K'^TP , the re- 
action mechanism for K'^YP is most certainly more 
complicated compared to K^K, and thus more dif- 
ficult to model correctly. 

2. The models reproduce reasonably well the forward 
peaking strength in ctj/, ctlt, and gtt for K^K 
and i^+E° for both final states for higher W . At 
W < 2 GeV where the resonance contributions are 
a larger contribution relative to the non-resonant 
background, the agreement is noticeably worse. 

3. None of the models reproduces the trends in ctlt' 
for either final state across the full W spectrum. 
Interestingly, the hadrodynamic model of Maxwell 
et al. that includes the available tlt' data from 
Ref. [13 has by far the worst agreement with these 
data, although the available g^t' data only go up 
to Q2=i.o GeV^. 

4. The GLV Regge model that includes no s-channel 
resonance terms, does as well as any of the other 
models in describing these data. For the K^YP 
final state for < 2 GeV, which has strong s- 
channel contributions, the GLV model significantly 
underpredicts <tu- However, for K, which has a 
much more significant i-channel exchange compo- 
nent within the resonance region, the GLV model 
underpredicts ajj for < 1.9 GeV. But for W > 
2.2 GeV, the GLV model well matches the data 
for both final states over our full kinematic phase 
space. 

5. For A'+A, the RPR-2011 model fares noticeably 
worse than for the RPR-2007 model over all angles 
for W < 2.1 GeV for all of the structure functions. 
For higher W, where the response is essentially 
fully i-channel, the RPR-2007 and RPR-2011 mod- 
els agree well with the data and with each other. 

B. Energy Dependence 

To more directly look for s-channel resonance evidence, 
the extracted structure functions are presented as a func- 
tion of the center-of-mass energy W for our ten values of 
cos^Ij-. Figs, [in] and [T7] show the results for our K^A 
and K+T,° data, respectively, at Q^=1.80 GeV^. 

A number of observations can be made regarding the 
data: 

1. For K^A production, <tu shows a broad peak at 
about 1.7 GeV at forward angles, and two peaks 
separated by a dip at about 1.75 GeV for our 
two backward angle points. This corroborates 
similar features seen in recent pho to- and elec- 
troproduction results [H, [56l458l |. Within ex- 
isting hadrodynamic models, the structure just 
above the threshold region is typically accounted 
for by the known 7V(1650)l/2-, A^(1710)l/2+, and 



iV(1720)3/2+ nucleon resonances. However, there 
is no consensus as to the origin of the bump feature 
at '--^1.9 GeV that was first seen in the if+A photo- 
production data from SAPHIR [s^. It is tempting 
to speculate that this is evidence for a previously 
"missing", negative-parity J = 3/2 resonance at 
1.96 GeV predicted in the quark model of Capstick 
and Roberts [4j] . This explanation was put forward 
in the work of Bennhold and Mart [59|, in which 
they postulated the existence of a 3/2~ state at 
1.9 GeV. However, in Ref. [2l| it was shown that a 
A^(1900)3/2+ state is required to explain the beam- 
recoil polarization data for K~^A. In Ref. [gO] this 
broad bump in the if+A cross section could be ex- 
plained by accounting for w-channel hyperon ex- 
changes. 

2. For K~^A, olt has about 20% of the strength of 
au and is consistently negative. For X+E°, glt is 
nearly zero everywhere except for VF=1.9 GeV at 
back angles. 

3. The (Jtt structure function is quite similar for 
K^A and K^TP over all kinematics with a strength 
comparable to olt- 

4. For JC+A, (Tlt' shows significant structure for W 
below 2.2 GeV. For higher W it is consistent with 
zero. 

5. In the K'^YP channel, ajj is peaked at about 
1.9 GeV, which also matches the photoproduction 
result [43, 113, HI] . CTTT, while small, shows a broad 
feature in this same region. These features are con- 
sistent with a predominantly s-channel production 
mechanism. In this region, beyond the specific N* 
resonances believed to contribute to AT+A produc- 
tion (and hence are strong candidates to contribute 
to K~^TP production), there are a number of known 
A* resonances near 1.9 GeV [l^ that can con- 
tribute to the K^YP final state, particularly the 
A(1900)l/2- and A(1910)l/2+. These A* states 
are forbidden to couple to the K^A state due to 
isospin conservation. 

The comparisons of the model calculations to the data 
clearly indicate that significant new constraints on the 
model parameters will be brought about when these new 
electroproduction data are included in the fits. We con- 
clude that the W dependence of A and K^YP produc- 
tion provides strong evidence for baryon resonance activ- 
ity within the reaction mechanism, but that the data in 
comparison to present models do not allow any simple 
statement to be made. We further conclude that at the 
current time the models that are limited to fits of the 
photoproduction data only, cannot adequately describe 
the electroproduction data. 
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FIG. 16: (Color online) Structure functions au, ult, (^tt, and ai^xi (in nb/sr) for K^A production vs. W at 5.499 GeV for Q^=1.80 GeV^ and for the 10 cosO}^ 
values. The error bars represent the statistical uncertainties only. The curves are defined in the caption of Fig. 1121 




FIG. 17: (Color online) Structure functions at/, olt, cttt, and oj^t' (in nb/sr) for K^Yf production vs. W at 5.499 GeV for (3^=1.80 GeV^ and for the 10 cos6*j^ 
values. The error bars represent the statistical uncertainties only. The curves are defined in the caption of Fig. 1121 
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C. Dependence 

Our data set provides a large reach and it is in- 
structive to study the W spectra for increasing values of 

. These data are shown in Figs. [TH] and [12] for the 
iiT+A and K^YP final states at two representative W 
points, 1.725 and 1.925 GeV. Included on these plots arc 
the photoproduction difi^erential cross sections for X+A 
from Ref. [11] and i<:+S° from Ref. at Q2=o for the 
kinematic points where they are iiyailable. Also shown 
are the data from gu from Ref. |16] from two different 
data sets, (i). Eb^2.h%l GeV, Q2^0.65, 1.0 GeV^ and 
(ii). £;b=4.056 GeV, Q^^i.q, 1.55, 2.05, 2.55 GeV^ at 
kinematic points that are reasonably close to the present 
data. 

What is seen by studying the evolution of ou is 
a reasonably smooth fall-off from the photon point. As 
the photoproduction data involve a purely transverse re- 
sponse, this smooth fall-off to finite in these kine- 
matics predominantly indicates a small longitudinal re- 
sponse. This is also indicated by the small strengths of 
OhT and ult' relative to ajj in Figs. [12] to [17] for back- 
and mid-range angles for the X+A final state and for all 
angles for the K^YF' final state. However, there is clearly 
a non-negligible longitudinal response in the iiT+A data 
at forward angles and for higher W as seen in these data 
(and also seen in the data of Ref. H^). Note that the 
comparisons shown in Figs. [TS] and [T^] are only for qual- 
itative comparisons as the kinematics are not a perfect 
match in all cases from Refs. [iB, EE, to the present 
data. 

The smooth fall-off of uu with increasing is consis- 
tent with the findings of the lower analysis of i^+A 
and K^YP electroproduction from Ref. [l^. As was the 
case in that work, it is seen that the interference struc- 
ture functions ctlt, cttt, and ctlt' for both final states 
do not demonstrate any strong dependence. How- 
ever, detailed comparisons with available models will be 
important to gain insight into the associated form factors 
for the iV* resonances found from fits to the photopro- 
duction data. 



D. Legendre Fits 

In order to investigate the possible evidence for the 
presence of s-channel resonance contributions in the sep- 
arated structure functions, we have considered two dif- 
ferent approaches. The first is with a fit of the individual 
structure functions ctc/, ctlt, <^tTi and ctlt' versus co^Q*^ 
for each and W point for the if+A and X+E" final 
states using a truncated series of Legendre polynomials 
as: 



U,LT,TT,LT' 



Pi{cos9K) dcose*j^. (21) 



The structures seen in these coefficients versus W are 
likely indicative of s-channel contributions. Note that 
the appearance of a structure at a given value of W in 
each of the different coefficients most likely suggests 
the presence of a dynamical effect rather than the signa- 
ture of an N* contribution. Instead, the appearance of a 
structure in a single d coefficient at the same W value 
and in each of the points is more likely a signal of an 
N* contribution. 

The fits for K+A show structures at GeV in 

Co for both au and a^x, W=1.9 GeV in C2 and C3 for 
au, and Vl^=2.2 GeV in C3 for cru. The fits for K+Y° 
show structures at Vl^=1.9 GeV in Co and C2 for au 
and (TTT- Of course, making statements regarding the 
possible orbital angular momentum of the associated s- 
channel resonances requires care as interference effects 
among the different partial waves can cause strength for 
a given orbital angular momentum value to be spread 
over multiple Legendre coefficients. 

In a second approach, each of the Legendre coefficients 
can be further expanded in terms of products of pairs 
of multipole amplitudes, but these expansions quickly 
become unwieldy as the number of participating partial 
waves increases. However, one simple thing that can be 
done for additional insight is to fit the structure functions 
with a coherent Legendre series of the form: 



dau, 



LT,TT,LT' 

~dn* 



Y,Ci{Q\W)P,{cose*„) 



c^ 



The fit coefficients for i' = 3 are shown for K 
in Fig. [10] and for K+YP in Fig. [H] for = 1.80 GeV^. 



(22) 

Here the Pi are the usual Legendre polynomials. The co- 
efficients Ci{Q^,W) are the amplitudes of the coherent 
5, P, and I?-wave contributions, respectively, while 
takes into account a incoherent "background" connected 
with higher-order terms that are not taken into account 
in the truncated sum. Of course, one must take care 
against making too much of the fit results using the sim- 
plistic form of Eq. (|22|) . This approach is not meant to be 
an attempt at a true amplitude fit. Rather the point is 
to look for structures that appear at a given W and for 
each for a given Ci coefficient as suggestive evidence 
for possible TV* resonance contributions. Fig. [H] shows 
the Legendre coefficient from this approach for au for 
the K^K reaction for the three points in this analy- 
sis. Fig. [53]is the corresponding figure for K^Y'^ . 

The fit coefficients for au shown in Figs. [H] and 
[23] show reasonable correspondence among all three 
points. For the if+A fits, strength is seen at: 
M^=1.7 GeV in Co, M^=1.9 GeV in Ci, and 1^=2.2 GeV 
in C2 . While it might be tempting to view this as corrob- 
oration of the findings of the if+A photoproduction am- 
plitude analysis from Ref. [6l|, obviously more detailed 
work is required. For the K'^Y^ fits, strength is seen 
at W=im GeV in Co and 14^=1.9 GeV in C2. It is 
interesting that there is no signature of strength in the 
P-wave as seen through the coefficient Ci, but again a 
higher-order analysis will be required to make more def- 
inite statements. 
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(K A-W=1.725 GeV) 



(K A-W=1.925 GeV) 



FIG. 18: (Color online) Structure function au vs. (GeV^) for the K+ k final state for two values of VF=1.725 and 1.925 GeV. 
The black circles are the data from this work, the red squares are the photoproduction points from Ref. 46] , and the green stars 
and triangles are from the lower data from Ref. [10I- The error bars include both statistical and systematic uncertainties. 
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FIG. 19: (Color onhne) Structure function au vs. (GeV^) for the K+Yp final state for two values of VK=1.725 and 
1.925 GeV. The black circles are the data from this work, the red squares are the photoproduction points from Ref. f47|, and 
the green stars and triangles are from the lower data from Ref. fig] . The error bars include both statistical and systematic 
uncertainties. 



VII. SUMMARY AND CONCLUSIONS 



We have measured if+A and K'^YP electroproduc- 
tion off the proton over a wide range of kinematics in 
the nucleon resonance region. We have presented data 
for the differential cross sections and separated structure 



functions ay, ctlt, frr, and ult' for from 1.4 to 
3.9 GeV^, W from threshold to 2.6 GeV, and spanning 
nearly the full center-of-mass angular range for the . 
In addition to the increased kinematic reach of these data 
relative to the previously published K^Y electroproduc- 
tion structure functions from CLAS in Ref. 16:], this new 
data set is an order of magnitude larger, allowing for finer 
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FIG. 20: Legendre polynomial fit coefficients (nb) from 
Ea. (|2ip vs. W for the K'^A separated structure functions 
ffu, cfLT, o-TT, and a^x' for (3^=1.80 GeV^. 
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FIG. 21: Legendre polynomial fit coefficients (nb) from 
Ea. (|21|l vs. W for the K^YP separated structure functions 
o"(7, (tlt, o-tt, and ai^-r' for — GeV^. 
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FIG. 22: Coherent Legendre polynomial fit coefficients 
((nb/sr)^/^) from Eq.dm vs. W for the K+A separated struc- 
ture function au for Q^^l.SO, 2.60, and 3.45 GeV^ 
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FIG. 23: Coherent Legendre polynomial fit coefficients 
((nb/sr)^/^) from Eq.^ vs. W for the K+T," separated 
structure function au for Q^=L80, 2.60, and 3.45 GeV^. 



binning in W and cos^^. 

The structure function data for both K^A and if+S" 
indicates that for W below 2.2 GeV and back angles, 
there is considerable strength of contributing s-channel 
resonances for K^A and if+E°. For higher W, the t- 
channel non-resonant background dominates and the re- 
action dynamics are well described solely through inter- 
ference of K and K* Regge trajectories. 

A Legendre analysis confirms these qualitative state- 
ments. For the K~^A final state, the Legendre moments 
of the structure functions indicate possible s-channel res- 



onant contributions in the 5- wave near 1.7 GeV, in the 
P-wave near 1.9 GeV, and in the D-w&ve near 2.2 GeV. 
This is in qualitative agreement with the more detailed 
amplitude analysis of Ref. [6l|. For the K^T,'^ final state, 
strong 5*- wave strength is seen at 1.8 GeV and strong D- 
wave strength is seen above 1.9 GeV, precisely where sev- 
eral A* states are expected to couple. Of course more de- 
tailed and quantitative statements await including these 
data into the coupled-channel partial wave fits. Such 
analyses would help to provide important complementary 
cross checks to the fit results of the recent Bonn-Gatchina 
coupled-channels results from Ref. 7] that seem to favor 
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a much richer mix of states to describe the available pho- 
toproduction data. 

Finally, detailed comparisons of our data have been 
made with several existing models. These include the 
hadrodynamic model of Maxwell et al. [13] that has been 
constrained by both the CLAS photo- and electroproduc- 
tion data sets (both cross sections and spin observables) , 
the Regge model of Guidal et al. 35] that has only been 
constrained by high-energy photoproduction data to fix 
the parameters of the Regge trajectories, and the Regge 
plus resonance model from Ghent [2^ that has been con- 
strained by the existing high statistics photoproduction 
data. None of the available models does a satisfactory job 
of describing the structure functions below W — 2 GeV 
for either K^A or if+E". In fact, several of the more 
recent models (e.g. RPR-2011 and the MX model includ- 
ing the CLAS ctlt' data) actually are in worse agreement 
with the data below 2 GeV than for earlier versions of the 
models. Clearly more work on the modeling and possibly 
the fitting/convergence algorithms is required to be able 
to fully understand the contributing N* — >■ ivT+A and 
iV*, A* — > if+E" states and to reconcile the results from 



the single-channels models with the currently available 
coupled-channel models. 
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